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ABSTRACT 

Determining fundamental properties of stars through stellar modeling has improved substantially 
due to recent advances in asteroseismology. Thanks to the unprecedented data quality obtained by 
space missions, particularly CoRoT and Kepler, invaluable information is extracted from the high- 
precision stellar oscillation frequencies, which provide very strong constraints on possible stellar models 
for a given set of classical observations. In this work, we have characterized two relatively faint stars, 
KIC 10920273 and KIC 11395018, using oscillation data from Kepler photometry and atmospheric 
constraints from ground-based spectroscopy. Both stars have very similar atmospheric properties; 
however, using the individual frequencies extracted from the Kepler data, we have determined quite 
distinct global properties, with increased precision compared to that of earlier results. We found 
that both stars have left the main sequence and characterized them as follows: KIC 10920273 is a 
one-solar-mass star (M = 1.00 ± 0.04M Q ), but much older than our Sun (t = 7.12 ± 0.47 Gyr), while 
KIC 11395018 is significantly more massive than the Sun (M = 1.27 ± 0.04M Q ) with an age close to 
that of the Sun (r = 4.57 ± 0.23 Gyr). We confirm that the high lithium abundance reported for 
these stars should not be considered to represent young ages, as we precisely determined them to be 
evolved subgiants. We discuss the use of surface lithium abundance, rotation and activity relations as 
potential age diagnostics. 

Subject headings: asteroseismology — stars: evolution — stars: fundamental parameters — stars: 
individual (KIC 10920273, KIC 11395018) — stars: solar-type 



1. INTRODUCTION 

Classical modeling of single stars mostly relies on fit- 
ting to the atmospheric properties obtained through 
spectroscopic and/or photometric observations, such as 
effective temperature, surface gravity, and elemental 
abundances. This yields a large number of possible mod- 
els covering a wide range of values for the fundamental 
properties, particularly in the absence of independent 
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radius and luminosity measurements. Asteroseismology 
has been revolutionizing stellar modeling as a result of 
the precise and accurate inferences on stellar structure 
that have been made possible by a new generation of as- 
teroseismic observations. Stellar fundamental properties, 
particularly mass and radius, can be determined to a few 
percent uncertainty, even when only the average seismic 
parameters are used as additional constraints, with the 
precision on age determination being as good as ^20%. 
The precision and accuracy of these properties increase 
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further when t he individual oscill ation frequencies are 
used (see, e.g., iMathur et al.l 120121 for a comparison be- 
tween using average and individual seismic quantities as 
constraints in the modeling of a large sample of Kepler 
stars). Moreover, the individual frequencies allow us to 
obtain information about the stellar interiors. 

Asteroseismology has proved very effective in con- 
straining the stellar age and the evolutionary stage with 
the help of specific features seen in the oscillation spectra. 
Mixed modes are particularly important in this regard. 
As a star evolves, the frequencies of the p modes decrease 
due to increasing stellar size, while the g-mode frequen- 
cies increase. By the time the star moves off the main 
sequence, the g- and p-mode trapping cavities are closer 
to each other, which results in the interaction of the two 
types of modes as they go through "avoided crossings" . 
The modes affected by this interaction are referred to 
as mixed modes due to having g-mode characteristics in 
the deep interior, and p-mode cha ract eristics near the 
surfa ce, of the star Csee lOsak"Hll975l and lAizenman et all 
119771 for an introductory discussion). They are sensi- 
tive to the central conditions and hence encode infor- 
mation from the core, where the chemical composition 
changes due to the nuclear reactions driving the evolu- 
tion of the star. Since the timescales of avoided cross- 
ings are very small compared to the stellar evolutionary 
timescale, mixed modes prov ide very strong constraints 
on the stellar age (see, e.g.. IDeheuvels fc Michel! 120111 : 
iMetcalfe et al.ll2010 ; Bcnom ar et al.ll2012l for recent anal- 
yses). 

Kepler is a space telescope with a diameter of 0.95 m 
that has been providing high-quality photometric data 
sinc e the beginning of its operation s in May 2009 (see, 
e.g..lBorucki et al.ll20Tfl iKoch et alJl2010t iChaplin et alJ 
1201011201 If ). The mission's primary objective is to search 
for Earth-sized planets through the transit method. 
Asteroseismology is being used to characterize a sub- 
sample of stars, some of which host planets. Kepler 
monitors more than 150,000 stars, and ^2000 of these 
were selected to be monitored for one month each 
in short-cadence mode (58.9 s int egrations) during the 
first ~10 months of the mission (IGilliland et al.l [20101 : 
IChaplin et al.1 120111 ). Solar-like oscillat ions were de- 
tecte d in at least 500 of those survey stars (|Chaplin et al.l 
1201 1|) . A subsample (~190) of these have been mon- 
itored for more than 3 months, and precise determi- 
nation of the oscillation properties has been com pleted 
for part of the sample (|Appourchaux et al.1 12012). As- 
teroseismology has been proving successful in determin- 
ing their global properties and inferring their interiors 
(see, e.g. IMetcalfe et al.|[2010l [20121: iCreeyev et all 12011 
IDeheuvels et alJl2012t IMathur et alJl2012t ). 

KIC 10920273 (kepmag=l 1.93 mag, i.e., apparent mag- 
nitude as observed through the Kepler bandpass) and 
KIC 11395018 (kepmag=l 0.76 mag) are among a hand- 
ful of asteroseismic targets that were observed continu- 
ously from the start of science operations. Consequently, 
extended timeseries were available from early in the 
mission, making b oth stars attractive targets for aster- 
oseism ic analysis (|Campante et al.1 120111 : IMathur et al.1 
1201 If ). We also acquired ground-based spectra in order 
to characterize these stars. They are G-type stars with 
very similar spectroscopic properties, especially T e g and 
logg (see, Sect. 2.1), so it is difficult to discriminate be- 



tween models for the two stars using classical approaches. 
These approaches include matching the position of the 
star in the Hertzsprung-Russell (H-R) diagram in the 
form that shows luminosity versus effective temperature 
as the star evolves, or alternatively in the \ogg-T e ff di- 
agram, given that the luminosity cannot be calculated 
using the available observations. 

We present the observational data employed to char- 
acterize our stars in Section 2, our modeling approach in 
Section 3, and the results in Section 4, while Section 5 
provides a summary and conclusions. 

2. OBSERVATIONAL CONSTRAINTS 
2.1. Atmospheric properties 

Atmospheric properties of KIC 10920273 and 
KIC 11395018 were obtained from observations w ith the 
FIES spectrograph (|Frandsen fc Lindbergl 119991 ) at the 
Nordic Optical Telescope (NOT on La Palma, Spain), 
at medium resolution (R~46,000) in July and August 
2010. The reduced spectra were analy zed by several 
teams and the results were presented by iCreevev et al.1 
(2012). The constraints we used for our analysis are 
shown in Table[T] (see also the 1- and 2-tr error boxes in 
Fig.[T]). We adopted a set of atmospheric constraints for 
each star that were closest t o the mean of result s from 
several methods described bv lCreevev et a l. (2012). This 
approach was preferred for the sake of reproducibility, 
rather than using the mean values. However, we did 
not restrict our model-searching space to less than 3-cr 
uncertainty around these constraints; therefore, the 
selected values represent well the overall results of the 
spectroscopic analysis. 

The spectroscopic usini values for KIC 10920273 and 
KIC 1139501 8 are 1.5 ± 2.2 km s" 1 and 1.1 ± 0.8 kms" 1 , 
respectively (|Creevev et al.ll201~2l ). These low values in- 
dicate either slow rotation or low inclination angle i, al- 
though the latter is statistically unlikely. Rotational pe- 
riods from modulation of the Kepler data due to spots 
on the surfaces of the two stars were measured to be 
~27 days for KIC 10920273 (iCampante et alJl20Tl . and 
~36 days for KIC 11395018 iMathur et alJl201lfc The 
signal-to-noise ratios (SNR) of the peaks used for these 
measurements, particularly for KIC 10920 273, were low, 
so the results should be used with caution. IMathur et al.l 
(1201 lh inferred i ■ 45° for KIC 11395018 based on com- 
bining the rotational frequency splittings with the mea- 
sured rotation period, which implies this star to be a 
slow rotator. A wider range of inclinations w as pos- 
sible for KIC 10920273 (|Campante et al.|[20TTI ). How- 
ever, the modulation in the light curve has been detected 
with less uncertainty in the new analysis performed using 
longer time series including Kepler Q9 and Q10 data of 
KIC 10920273 (Garcia et al., private communication). If 
confirmed, this would be consistent with a relatively high 
inclination angle, and slow rotation. The effects of the 
centrifugal force are negligible for slowly rotating stars. 1 
However, rotational mixing may lead to changes in the 

1 If, counter to our expectations, one of these stars were to be 
confirmed as a fast rotator, rotational effects on t he oscillation 
frequ encies would have to be taken into account (see lSuarez et at] 
12010 for a detailed analysis of the effects of centrifugal distortion 
on solar-like oscillations). Currently, there is no robust detection of 
rotational frequency splittings that can be included as constraints 
in our modeling. 
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TABLE 1 

Adopted atmospheric cons traints for KIC 1092 0273 
and KIC 11395018 (from ICreevey et al1|20T2T) 



Star 


Tcff(K) 


logs 


[Fe/H] 


KIC 10920273 


5790 ± 74 


4.10 ±0.10 


-0.04 ±0.10 


KIC 11395018 


5700 ± 100 


4.10 ±0.20 


0.13 ±0.10 



properties of the models even for slowly rotating stars, 
because the efficiency of this mixing is more directly re- 
lated to differential rotation in stel lar interiors rather 
than to surface rotational velo cities (Pinsonn eault et"al"1 
Il990t lEggenberger et al.ll2010f) . Studying the impact of 
rotation on post-main-sequence stars would require a de- 
tailed discussion of the effects of rotational mixing on 
the chemical gradients in the central parts of the star. 
These influence the asteroseismic properties of the mod- 
els, particularly the mixed modes. In the specific case of 
the evolved post-main-sequence stars modeled here, how- 
ever, we expect that these effects on the ch emical gradi- 
ents w ould already be erased, as found bv iMiglio et al.l 
(|2007h for models of 12 Bootis A in the thick-shell- 
H-burning phase (see also the discussion of the effects 
of microscopic diffusion i n the subgiant HD 49385 by 
iDeheuvels &: Michell 1201 ID . We therefore have not in- 
cluded the rotational effects for most of the analyses (see 
Section 3). 

When the atmospheric properties alone are considered, 
these two stars are very similar. Due to the degen- 
eracy inherent in the H-R diagram analysis (see, e.g., 
Fernandes fc Monteiroll2003() . it is not possible to deter- 
mine the global stellar properties with sufficiently high 
precision to study their detailed characteristics without 
the help of seismic data, which we now discuss. 
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Fig. 1.— Logg-T cff diagram for KIC 10920273 and 
KIC 11395018. Surface gravity, g, is in cgs units. Spectroscopic 
constraints given in Table [T] are shown by 1- and 2-ct error boxes 
(dotted; blue for KIC 10920273 and red for KIC 11395018). Evolu- 
tionary tracks of two models indicated by star symbols (SA1 and 
BA1 from Tables[4] and [5j are plotted using the same color code. 
The points with error bars represent the weighted means and the 
standard deviations of the asteroseismic determinations (see Ta- 
bles[4]and[5]l. 



2.2. Asteroseismic data 

We used Kepler data from observations made in the 
period from 2009 May to 2010 March, i.e., from the com- 
missioning run (Q0) through Quarter 4 (Q4). The formal 
freq uency resolution i s ~0.05 jnHz. From the power spec - 
tra, iCampante et all (|2011[ ) and iMathur et all (|201lD 
reported individual frequencies for KIC 10920273 and 
KIC 11395018, based on analyses performed by several 
teams. The final sets of results included a minimal and a 
maximal list of frequencies, where the former were those 
agreed upon by more than half of the fitters and the latter 
were those agreed upon by at least two fitters. There- 
fore, the frequencies that are in the maximal list but 
not the minimal list are less certain. For details of the 
frequency-extraction techni ques and the sele c tion m eth- 
ods, we refer t he re ader to ICampante et al.l (|201lD and 
IMathur et al.l (|2011[) . The analysis of each star resulted 
in the extraction of up to a total of 25 individual os- 
cillation frequencies for radial (/ = 0), dipole (I = 1), 
and quadrupole (I — 2) modes, including several mixed 
modes. These mixed modes carry information from the 
core and hence provide stronger constraints on the evolu- 
tionary stage of the stars, as discussed in Section 1. We 
started by searching for models using the minimal-list 
of frequencies and then extended our analysis to include 
additional frequencies from the maximal lists. 

3. MODELING APPROACH 

Asteroseismic modeling is performed by optimizing the 
stellar model parameters to match the observed seismic 
quantities and also the classically observed (or derived) 
stellar properties, such as effective temperature, surface 
gravity, surface metallicity (along with radius, mass and 
luminosity, when available). The seismic quantities in- 
clude, but are not limited to, the average large and small 
frequency separations 2 , the frequency of maximum power 
in the oscillation spectrum (f ma x), and the individual os- 
cillation frequencies. Naturally, the individual frequen- 
cies provide the most detailed information and the high- 
est precision in the deri ved stellar properti es (see, e.g., 
iMetcalfe et~aT1l2010l and IMathur et alJl2012t ). 

We used the individual oscillation frequencies and the 
atmospheric properties (T c g, logg, and [Fe/H]) as con- 
straints to carry out the stellar model optimization. As 
an initial guess for the parameter space to be searched, 
we used preliminar y results of th e mass determination 
from the analysis of lCreevey et al.l (j2012h , which were in 
agreement with their final results within the uncertain- 
ties. They derived stellar properties using the average 
seismic quantities together with the atmospher i c con - 
straints. The final values given bv lCreevevetaLl ([2012T ) 
were 1.25±0.13 M Q for KIC 10920273 and 1.37±0.11 M e 
for KIC 11395018. 

Five teams participated in the modeling of these stars 
using a variety of evolutionary codes and fitting methods. 
Most of the methods were cither based on searching for 
the best-fitting model in a grid specifically computed for 
this analysis or on using a pre-existing grid to determine 
the general area of the stellar properties in the parame- 

2 The large frequency separation is Au n ; = u n ; — ^ n _i,i, and 
the small frequency separation is Sv n ; = v n i — u„_i ;_|_2, where 
v n i is the frequency of the mode with spherical degree I and radial 
order n. 
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ter space before going into further refinement process for 
individual stars. One team used the Asteroseismic Mod- 
eling Portal (AMP), which is a pipeline analysis tool that 
optimizes the seismic and non-seismic properties globally 
using a genetic algorithm. AMP starts the model-search 
with four random independent sets of initial parameters 
and performs the se a rch over a large para meter space 
([Metcalfe et al.|[200l IWoitaszek et aij|2009Q . The vari- 
ety of codes and methods employed give us an estimate of 
the external uncertainties inherent in the analysis. The 
list of codes and the configurations regarding the input 
physics are presented in Table [2] 

The individual fitting methods also differed slightly. 
ASTEC1 calculated grids of models within the 3-er un- 
certainty of the non-seismic constraints and performed 
the optimization by a 2-step process, refining the grids 
several times in the second step guided by the seismic 
X 2 values - described by Equation (J3J. ASTEC2 ex- 
plored the models, which included turbulent diffusion, 
and calculated individual models guided by the frequen- 
cies. CESAM looked for models reproducing the first 
avoided crossing as an initial requirement and then per- 
formed an optimization using y -minimization to deter- 
mine stellar mass and age (see iDeheuvels fc Michell[2011l 
for details of this method) . The Geneva stellar evolution 
code was used to compute grids of rotating models with 
an initial velocity of 50kms _1 on the zero-age main se- 
quence (ZAMS). This value results in surface velocities 
that are typically lower than lOkms -1 at the end of the 
main sequence (MS) for a solar-type star that is assumed 
to undergo magnetic breaking on the MS due to the 
prese nce of a convective envelope ([Krishnamurthi et al.1 
1997). The initial parameters used by each team are 
given in Table [3J 

Oscillation fre quencies of low-degree modes were calcu- 
lated by LOSC ([Scuflaire et all 120081) for stellar models 
computed by CESAM, w hile the Aarhus Adiabatic Pul- 
sation Package (ADIPLS, Christcnscn-Dalsgaard 2008b) 
was used to calculate the frequencies for all of the other 
models. 

In relation to the oscillation frequencies, there is a 
well-known offset between the observed and the model 
frequencies for the Sun and solar-type stars, due to in- 
accurate representation of the near-surface layers in the 
models. To address t his issue, we us e d the empirical cor- 
rection suggested bv lKieldsen et al.l (|2008f) . who showed 
that the difference between the observed and calculated 
solar frequencies, which gets larger with increasing fre- 
quency, can be fitted by a power law: 



^obs(«, 0) - ^bost("-, 0) = a 



( i/ obs (n,0) 



V() 



(1) 



where v ^(ri, 0) and fb e st(n> 0) are the observed and best 
model frequencies with degree I = and radial order 
n, Vq is a constant frequency usually chosen to be the 
frequency at maximum oscillation power, a is the size of 
the correction at v and can be calculated for each model, 
and b is the exponent to be determined. 

The right- hand- side of Equation ([T]) is the correction 
term to be added to the acoustic (p-mode) frequencies 
of the best models. The mixed modes, however, are less 
sensitive to the properties of the near-surface layers since 
much of their energy is confined to the stellar center. In 



other words, we need to apply a smaller near-surface cor- 
rection to the frequency of a mixed mo de than to a p 
mode with a similar frequency. Following Brandao et al. 
(i2on . we scaled the magnitude of the correction in- 
versely with Q n i , the inertia of a given mode normalized 
by the inertia of a radial (I = 0) mode at the same fre- 
quency (see, e.g. JAerts et al.1120101 ). Note that the inertia 
of a mixed mode is much higher than that of a p mode. 
The correction to be applied to all calculated frequencies 
is then of the form: 



*Wr(n, = fcw(rc, l) + a 



Q n 



(2) 



where ^ CO rr represents the corrected model frequencies. 
Notice that v o \> s (n,0i) on the right-hand-side of Equa- 
tion ([T]) is replaced by the best model frequencies in order 
to allow us to correct the freq uencies outside the ra nge 
of observed radial modes (see, IBrandao et al.l [2011L for 
details) . 

The solar value of th e exponent b was calculated 
by IKieldsen et all <f200l to be 4.90 using the GOLF 
data (lLazrek et al.ll 199 71) and the solar "Model S" of 
IChristense n-Dalsg aard et al.l (fl9 96) . This value was 
found to range from 4.40 to 5.25 for the same model, 
depending on the number of radial orders included in 
the calibration, but a was found to v ary b y less than 
0.1 ^Hz in all cases. IKieldsen et al.l ([2008D suggested 
that the solar b value may be used for solar-like stars 
and this approach w as successfull y applied to Hyi 
(IBrandao et all I20TTI) HP 49385 (IDeheuyels fc Michel 
120111). KIC"TT02676 4 (IMetcalfe et all 120 101) and a sam- 
ple of Kepler stars (Mat hur et al.l 120121 ). In this work 
ASTEC1, ASTEC2, and Geneva codes ado pted the solar 
value b = 4.90 from IKieldsen et al.l ([2008D for calculat- 
ing the correction term, while AMP adopted b = 4.82, 
which is th e solar-calibrated valu e for AMP with the Bi- 
SON data (|Chaplin et al.|[l999l) . and CESAM adopted 
b = 4.25, the ca librated value using the GOLF data 
(|Gellv et al.|[200l . Given how little a varies for a rel- 
atively large range of b for a given model, as discussed 
above, using slightly different b values for model-fitting 
has a negligible impact on the results. 

4. RESULTS AND DISCUSSION 

4.1. Global properties 

To select the best models, we defined the two normal- 
ized x 2 measures shown in Equations ([3]) and Q, which 
allowed us to evaluate the qualities of the fits for the 
atmospheric parameters and the seismic parameters sep- 
arately. The seismic measure was 



Xscis 



N ^ 

n.l 



v bs(n,l) - v co „{n,l) 
a(v ohs (n,l)) 



(3) 



where N is the number of observed frequencies, Vcorrin, I) 
represents the near-surface-corrected model frequencies 
with spherical degree I and radial order n, Vobsin, are 
the observed frequencies, and cr(i/ bs(^, 0) are * ne uncer- 
tainties on the observed frequencies. The measure for 
the atmospheric properties was 



Xatm 



1 



obs 



c(Pobs) 



(4) 
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where P={T c g, \ogg, [Fe/H]} and the subscripts "obs" 
and "mod" represent the observed and model properties, 
respectively, with cr(P bs) denoting the observational un- 
certainties. The values of [Fe/H] for the models were 
calculated using the formula [Fe/H] = log(Z/A) mod — 
log(2?/X)f^, where the solar value was adopted from 
IGrevesse fc Noeli (fl993l) as (Z/X) e = 0.0245. 

Each modeling team returned the model that best 
matched the observational constraints. We present the 
properties of these models in Tables |4] and 03 along with 
the normalized % 2 values 3 . We also present models fitted 
using more or fewer frequencies than those in the mini- 
mal lists in order to see whether the model-fitting results 
change considerably. In each case, we calculated the Xg eis 
in the tables using only the frequencies that were com- 
mon constraints for all of the models, in order to achieve 
a consistent evaluation of the models. 

There is a good agreement between the observed fre- 
quencies and the model frequencies. The quality of the 
fits can be seen in the echelle diagrams for a sample of 
models (Fig. [2]). Overall frequency patterns, including 
the dipolar mixed modes, are matched quite well. The 
fact that Xseis>l-0 implies that either the observational 
uncertainties are underestimated or the models are in- 
complete representations of the observational data. The 
models with relatively high Xseis (~ 20.0), are those that 
either could be improved with further refinement or that 
do not reproduce all of the modes simultaneously, in par- 
ticular the mixed modes, which are more difficult to fit. 
Due to their strong sensitivity to stellar evolution, mixed 
modes tend to dominate the model-fitting results. The 
timescale on which the signatures of these modes evolve 
is very short, so it becomes difficult to find good fits un- 
less the grid of models used is very fine. 

We present our results in Tables |4] and [5] Properties of 
all the models are listed, along with the weighted mean 
values and the standard deviations. Both stars have 
left the main sequence (central hydrogen mass fraction 
X c = 0.0) but have quite different characteristics, as seen 
in Tables Hand El The fact that KIC 10920273 is an old 
solar analog (with one solar mass and near-solar metal- 
licity) makes it an interesting target for further studies. 
A typical solar model would turn off from the MS at 
around 9-10 Gyr. However, the metallicity and partic- 
ularly the helium abundance alter this age estimate. In 
this case it is the high helium abundance that affects 
the MS turn-off age more than the low metallicity. The 
models with higher helium abundance behave similar to 
those with higher mass (higher luminosity), following an 
evolutionary track similar to that of a higher-mass star, 
and hence have shorter MS lifetimes. 

The results in Tables H] and [5] were weighted by the 
goodness of the seismic fit, i.e. the inverse of Xs Gis . This 
way, any misleading contribution due to coarse grids is 
eliminated. Although this does not provide a direct mea- 
surement of the systematic uncertainties, it still allows us 
to estimate the order of magnitude of the external errors 
expected from using different inputs, codes, and fitting 
methods. A similar determination of the systematic er- 

3 The labels of the models presented in Tables 4 and 5 start with 
"S" and "B" for KIC 10920273 and KIC 11395018, which stand for 
"Scully" and "Boogie" — the nicknames of the stars within Kepler 
Asteroseismic Science Consortium, Working Group 1. 



rors for the case of bright Kepler s tars 1 6 Cyg A and B 
was carried out by Metcalfe et al.l (|2012[ ) using different 
evolutionary codes and fitting methods. The uncertain- 
ties we determined are mostly greater due to the lower 
SNR in the data of our faint stars. Systematic errors 
in determination of stellar properties using grid-based 
pipelines caused by different observational constraints 
and different input physics w ere discussed more g ener- 
ally for a few Kepler stars by iCreevev et al.l (|2012t ). We 
discuss the uncertainties further in the next section. 
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Fig. 2. — Echelle diagrams for the selected models of 
KIC 10920273 (SA3, upper panel) and KIC 11395018 (B P, lower 
panel ) . The radial orders are plotted horizontally (see IBeddingl 
2012) and the vertical axes indicate the frequencies at the middle of 
each order. The background is a gray-scale map of the power spec- 
trum from the observations, with the observed frequencies shown 
by filled symbols with l-cr error bars. Minimal-list frequencies are 
plotted in red, while the additional maximal-list frequencies are 
plotted in blue. Open symbols are used for the model frequencies, 
with smaller size implying larger normalized mode inertia. Circles, 
triangles, and squares represent I = 0, I = 1, and 1 = 2 modes 
respectively. 



4.2. Discussion of uncertainties 

To evaluate the typical uncertainties caused by dif- 
ferent input physics further, we calculated some addi- 
tional models and small grids using KIC 10920273 as a 
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test case. We selected SB1 as our base model. Keeping 
the input parameters (mass, Z/X, Y, a) fixed, we first 
explored the effects of changing one single ingredient of 
input physics at a time. We also changed the value of a 
while keeping everything else fixed. On every new evo- 
lutionary sequence, we calculated the frequencies of the 
models that had all atmospheric properties (log<?, T e s, 
and [Fe/H]) within 3 a of the observed values. We then 
selected the model that best matched the observed fre- 
quencies. In the cases of including core overshoot, chang- 
ing the convection treatment (to CG M formulation), and 
using severa l different nuclear rates (jBahcall et al.lll995t 
Parker 1986; Adclbcrger et al. 1998), the new models re- 
produced the observed frequencies as well as the base 
model with an age difference of only 0.7% at most (which 
corresponds to changes in radius of <0.2% and in log g of 
<0.1%). This means that the uncertainties on the final 
parameters caused by the corresponding inputs were neg- 
ligible. We note that the effect of core overshoot would 
be more significant for a main-sequence star. 

For the cases where we did not obtain a model having 
a seismic \ 2 comparable to that of the base model, we 
carried on with the analysis. These cases resulted from 
including diffusion and gravitational settling of helium, 
using two different versions of the low-temperature opac- 
ities given in Table [3J and varying the value of mixing 
length parameter (in the range of 1.6-2.0). For each of 
these cases, we computed additional small grids around 
the base model by varying all the input parameters, in 
order to see how much the output properties were differ- 
ent for two models with different input physics but sim- 
ilar frequencies. We then calculated the weighted mean 
and standard deviation in the same way as in the orig- 
inal analysis. The standard deviation in this case rep- 
resents the typical uncertainties caused by using a fixed 
set of input physics, hence decreasing the level of model- 
dependance in the results substantially. The mean val- 
ues for age, luminosity, radius, T c ff, and logg calculated 
from the additional analysis agreed with the original re- 
sults within 1 a, while their standard deviations were of 
the same order as those qiven in Table 2J This confirms 
that the uncertainties presented here are realistic. More- 
over, the resulting values of radius and logg from the 
additional analysis are essentially the same as the orig- 
inal results. This is reassuring given the importance of 
asteroseismology in determining the radius in a robust 
way. 

The internal uncertainties were different for each 
method. However, the dominant source of uncertainty is 
the non-uniqueness of the solution rather than the sta- 
tistical errors. Parameter correlations allow a trade-off 
between parameters, leading to different families of so- 
lutions that are almost equally good seismic fits. The 
effective range of these correlations is narrowed substan- 
tially, but not eliminated, by the use of seismic data. We 
used a single method (AMP) to evaluate the uniqueness 
of the best models for both stars, since the uniqueness 
depends on the specific constraints adopted in each case. 
AMP finds the lowest value of \ 2 in the entire search 
space. Along the way it also identifies the secondary 
minima - which can be far away from, but not much 
worse than, the formally best solution. For both stars, 
these secondary minima (SA2 and BA2 in Tables[4]and[5]) 
are marginally worse than the best solutions found using 



the same set of inputs (SA1 and BA1), but the values 
of the mass and helium abundance are quite different. 
This reflects the well-kno wn mass-He degeneracy (see , 
e.g.. ILebreton et al. l [l99l iFernandes fc Monteirol 120031: 
IMetcalfe et al.ll2009f) . which is not entirely lifted, even 
with the help of asteroseismology. The mass is con- 
strained more strongly than it would be without aster- 
oseismic data but in the absence of external constraints 
on the helium abundance, we cannot choose one partic- 
ular model. Consequently, we used both the primary- 
and secondary-minimum solutions from AMP in the cal- 
culation of the weighted means. The range of the AMP 
results leads to relatively large "uniqueness uncertain- 
ties" in the fitted stellar properties. These are deter- 
mined as follows for KIC 10920273, and KIC 11395018, 
respectively: 3%, and 7% in stellar mass; 8%, and 21% 
in the initial metallicity (Z/X)i\ 10%, and 25% in Y]; 
7%, and 8% in age; and 1%, and 2% in radius. Note that 
the larger uncertainties in the mass and chemical com- 
position for KIC 11395018 reflect the fact that it has two 
observed avoided crossings, which provide more strin- 
gent constraints such that neighboring models are signif- 
icantly worse. Only a relatively large jump along the pa- 
rameter correlations yielde d a secondary minimu m with 
a comparable seismic fit (cf. IMetcalfe etalll2010h . When 
we evaluate the uncorrelated statistical uncertainties for 
each of these minima with a local analysis using singu- 
lar value decomposition (SVD), we found very small er- 
rors (e.g., as low as 5xl0~ 3 % for mass), reflecting the 
steep x 2 surfaces corresponding to these results. Al- 
though these 'local' uncertainties are real, i.e., changing 
the values of the parameters by the correlated 'tiny' un- 
certainties causes a large increase in the x 2 (>50), the 
minima are not unique. Hence, a local analysis does not 
reflect the true uncertainties of the overall results. 

Inclusion of the maximal list of frequencies as con- 
straints in the analysis (five additional frequencies for 
KIC 10920273 and two for KIC 11395018) made the 
agreement between the model and observations better for 
both KIC 10920273 (SA3 in Table @J| and KIC 11395018 
(see BA3 in Table [5]). We also used a third list for 
KIC 10920273 to ensure that the less-certain frequen- 
cies did not bias our model-fitting. For that purpose, 
we left out two frequencies in the analysis. One of 
these was identified as a dipole-mode frequency belong- 
ing to the minimal list, but tagged as being close to 
the second harmonic of the long-cadence frequency (with 
v = 1135.36 ±0.31 /iHz), while the other one was identi- 
fied as a quadrupole mode (with v = 873.10 ± 0.32 fiRz) 
that was t agged as a possible mix ed mode introduced a 
posteriori (jCampante et al.l l201lT) . Both SA4 and SB2 
were selected using this alternative frequency set for 
KIC 10920273, and the agreement between the model and 
observations was not affected substantially. Therefore, 
we cannot ascertain whether these two peaks are stellar 
in origin. Nevertheless, we do not completely rule out 
the possibility of these peaks being stellar as some of the 
models that result from using alternative frequency sets 
do contribute to the weighted mean values significantly. 
Furthermore, we note that the peak at v — 873.10 /iHz, 
which is tagged as a possible quadrupole mixed mode, is 
in the middle of the frequencies of a dipole mixed mode 
and a quadrupole mode in most of our models. So the 
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observed peak may correspond to a dipole mixed mode, 
which, according to the models, has a relatively low nor- 
malized mode inertia indicating an observable amplitude. 

4.3. Comparison with previous results 

Comparing our results with those from the pipeline 
analyses of Creevey et al. (also given here in Tables 2] 
and [5j , we see that the mass determinations from the 
pipeline analyses were higher, which led to lower age es- 
timates. We emphasize that the previous pipeline analy- 
ses used only the average seismic quantities, hence lack- 
ing additional information from the individual frequen- 
cies and being affected by the uncertainties of the scal- 
ing relations. Therefore, some deviation from their val- 
ues was expected. Nonetheless, it is reass uring that the 
mass, radius, and age determinations of I Creevey et al.1 
(2012) are within 2-a uncertainty limits of our results 
for KIC 10920273, and within 1 a for KIC 11395018. We 
also confirm the robust determination of log g using scal- 
ing relations and grid-based analyses (see Table 7 in 
iCreevev et al.|[2012T) , with which our results are in agree- 
ment within 1 a. Additionally, we note that our results 
confirm that the mass and radius determined u sing only 
the scaling relations (e.g., iMathur et al.l I2012T ) provide 
good initial estimates for these properties (1.06±0.20Af© 
and 1.8O±O.11J? for KIC 10920273; 1.31±0.25M Q and 
2.21 ± 0.14i? Q for KIC 11395018). 

There is excellent agreement, for both stars, between 
our re sults and the mass estimates of Bcnomar et al.1 
(2012), who used the coupling strength of the observed 
mixed modes to determine the masses of several sub- 
giants, including KIC 10920273 (1.04 ± 0.04(±0.04)M Q ) 
and KIC 11395018 (1.21 ± 0.06(±0.04)M Q ). 

The most substantial improvement in this work comes 
from the use of individual frequencies which yield in- 
creased precision, with age being affected the most. The 
presence of the mixed modes in the data allowed us to 
determine the age with 5-7% precision, although with 
some model-dependency. This result is a major improve- 
ment on the 35-40% precision in age achieved using at- 
mospheric and mean seismic parameters. Both stars are 
determined to be post-main-sequence subgiants with no 
hydrogen left in their cores. Evolutionary tracks of the 
selected models are shown in Fig.[TJ 

Although we did not restrict the parameter search 
to be within l-a uncertainty around the spectroscopic 
constraints, the weighted mean values of T g from the 
models are within l-a limit for both stars, while log g 
results are in agreement with the spectroscopic values 
within 2 a (see Fig.QJ, and [Fe/H] within 1.5 a. Our 
log g results are in excellent agreement with asteroseis- 
mi c log q values obtain ed from scaling relations (given 
bv ICreevev et al.1 (|2012l ) and also in TablesH and [5]) . We 
also note that our temperature results are in good agree- 
ment with the revised photo metric values for the Keple r 
Input Catalog (KIC) from iPinsonneault et al.l (|2012j) . 
who derived T eS = 5872 ± 70 K for KIC 10920273, and 
T cff = 5650 ± 59 K for KIC 11395018. 

4.4. Non-seismic age diagnostics 

We discussed the asteroseismic constraints on the stel- 
lar age in Section 1. Here we discuss the implications of 
rotation and stellar activity on the age, as well as those 
of the surface lithium abundance. 



Rotation and activity are potentially valuable diagnos- 
tics of stellar age. It was shown that the Ca + emission 
luminosity, an indicator of stellar activity, dec ays roughly 
as t~ 0,5 for some cluster stars and the Sun (jSkumanichl 
11972t h furthermore, rotational decay was shown to follow 
the same law. Large samples of stellar rotation periods 
have been collected, and the Kepler mission promises 
many more. There is therefore substantial interest in 
stellar rotation-m ass-ag e , or gyrochronology, rel ations 
(see iBarnesI 120031 l2007t iMamaiek fe Hillenbrand! 120081: 
iMeibom et al.l 1201 lUEpstein fe Pinsonneaultll2012t ). We 
discussed in Section 2.1 that relatively slow rotation is 
inferred for both stars. Slow rotation rates imply rel- 
atively old stars, which is consistent with our astero- 
seismic determinations. However, one would not expect 
main-sequence spin-down relationships to apply directly 
to the evolved stars. Therefore, we cannot use the ro- 
tation rates for these stars to infer their ages with the 
age-rotation relations established for MS stars; these re- 
lations need to be calibrated for more evolved stars using 
larger samples. 

We have analyzed the chromospheric activity in the 
Can HK lines and found the levels of activity in both 
stars to be very low. Fig. [3] shows the Call K and H lines 
of KIC 10920273 and KIC 11395018 compared to the Sun. 
The solar spectrum was obtained from the solar light re- 
flected by Ganymede, which was observed with HARPS 
in April 2007 4 , when the Sun was close to the minimum 
of its activity cycle. We accounted for the different re- 
solving power of HARPS (R ~ 120,000) compared to 
FIES spectrograph (R ~ 46, 000) by convolving the so- 
lar spectrum with a Gaussian kernel of the appropriate 
width. It is clear that these two Kepler stars have chro- 
mospheric activity levels comparable to the quiet Sun, or 
lower. These low activity levels are consistent with the 
old ages we infer from asteroseismology; however, the 
rough nature of the empirical age- activity relations for 
post-MS stars does not allow us to make a quantitative 
analysis to infer ages. 

Another independent determination of stellar age may 
be obtained by measuring the Li content at the stel- 
lar surface. Lithium is easily destroyed in stellar inte- 
riors and is only produced under unusual circumstances; 
it has therefore been employed as an age indicator for 
low-mass stars. Lithium can be directly depleted if the 
surface convection zone is deep enough. It can also be 
mixed into the radiative interior, or it can be stored be- 
low the surface convection zone by microscopic diffusion 
processes. In standard stellar models, pre-main-sequence 
depletion occurs for most low-mass stars when they have 
deep conve ction zones, and it is most severe in lower 
mass stars (jBodenheimerl fl965[ ) . In a qualitative sense, 
a detection of Li in very cool stars is a strong indica- 
tor of youth. However, standard models also predict 
that stars of order 0.9 solar masses and higher would 
not experience main-sequence Li depletion, and there is 
strong evidence from open clusters for a steady decrease 
in Li as a functio n of time, even for stars more ma ssive 
than the Sun (s ec Zap palalll972l : lPinsonneaultlll997l : and 
ISestito fc Randichl 120051 for reviews. 

It was also shown bv lRandichl (|2010ft that, for a fraction 

4 http : //www. eso . org/sci/f acilities/lasilla/instruments/ 
harps/inst/monitoring/sun .html 




Fig. 3. — Chromospheric activity in the Call K and Call H lines (flux relative to the continuum) for KIC 10920273 (upper panels) and 
KIC 11395018 (lower panels). The solar spectrum (Ganymede taken in 2007 with HARPS) is overplotted with a dotted (red) line. The 
residuals between the stellar and solar spectra (at the bottom of each plot) show that the two stars have activity levels comparable to the 
quiet Sun or lower. 



of solar-like stars with effective temperatures between 
5750 K and 6050 K, Li is not further depleted after the 
age of ~1 Gyr, unlike for the Sun and many other stars. 
Due to this bimodal pattern, Li abundance alone cannot 
be used to determine the age for all stars, and a high Li 
abundance can only help define a lower limit for the age, 
since it may correspond either to a young star that has 
not yet depleted much Li, or to an older star that has 
st opped depleting Li a long time ago. 

ICreevev et al.l (|2012l ) showed that the two stars consid- 
ered here have strong Li absorption lines, which implies 
a high Li content at the surface (log N(Li)= 2.4 ±0.1 for 
KIC 10920273 and log N(Li)= 2.6±0.1 for KIC 11395018; 
where log N (Li) = log [n(Li)/n(H)] + 12, with n being the 
number density of atoms and log iV(H) = 12 by defini- 
tion). C onsidering the emp i rical L i-age relation estab- 
lished bv lSestito fc Randichl (|2005D . Creevey et al. then 
determined that the given Li abundances would indicate 
low ages (1-3 Gyr for KIC 10920273 and 0.1-0.4 Gyr for 
KIC 11395018), which are incompatible with the astero- 
seismic ages they determined through the pipeline mod- 
eling (see Tables 0] and [SJ performed using the average 
asteroseismic quantities. 

However, in addition to the bi-modality mentioned 
above, the age-Li relation has been shown to be valid 
for MS stars and does not necessarily extend to more 
evolved stars. This makes the age determination using 



the Li abundance ambiguous. Thus, despite the high 
Li abundance, we are confident that these are indeed 
evolved stars that have left the main sequence, due to the 
presence of the mixed modes in the observed oscillation 
spectra and as confirmed by our asteroseismic analysis. 

5. SUMMARY AND CONCLUSIONS 

We performed asteroseismic modeling of two Kepler 
stars, KIC 10920273 and KIC 11395018, for which we 
have long seismic data sets (>8 months) and ground- 
based follow-up spectra. We used individual oscilla- 
tion frequencies and atmospheric properties as initial 
constraints. We employed several evolutionary codes 
with different input physics, and various fitting meth- 
ods to determine the global stellar properties and esti- 
mate their uncertainties (see Tables 0] and [5]) . The near- 
surface correction was applied to the models, which re- 
produced the individual observed frequencies with con- 
siderable success; see Fig. [2] for a qualitative representa- 
tion with echelle diagrams. These two relatively faint 
stars, which have similar atmospheric properties accord- 
ing to the ground-based data, turned out to be substan- 
tially different - more than could have been predicted 
from their different metallicities - after incorporating the 
high-precision asteroseismic data into the modeling. 

KIC 10920273 resembles an old Sun, having one solar 
mass (1.00±0.04M Q ) and an age of r = 7.12 ±0.47 Gyr, 
while KIC 11395018 has a mass of 1.27±0.04M ra and an 
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age very close to that of the Sun (r = 4.57 ± 0.23 Gyr). 
These results agree, at the 2-cr level for KIC 10920273 
and 1-ct level for KIC 11395018, with the properties de- 
termined using the average asteroseismic quantities. The 
results presented here are much more precise than those 
from the average asteroseismic quantities, though, and 
they are also more accurate as a result of using more 
observational information, i.e. individual frequencies, 
and stronger constraints extracted from the observations, 
such as the mixed modes. We confirm the robust deter- 
mination of log g from the average seismic quantities, as 
our results are within 1-er uncertainty of the pipeline re- 
sults. 

We confirmed these stars to be subgiants (having 
evolved off the main sequence) and this allowed us to 
resolve the disagreement between the seismic ages deter- 
mined from the pipeline analyses and the ages estimated 
using the lithium abundance and the empirical Li-age 
relationship. Basically, the Li abundance cannot be em- 
ployed to estimate the ages of the subgiants. Similarly, 
existing age-rotation-activity relations can only be in- 
dicative for subgiants as these relations are calibrated 
mostly for the main-sequence stars. This must be taken 
into account for gyrochronology studies. 

We will soon obtain longer data sets from Kepler for 
many more stars and our results are a good indication of 
what we can achieve. We note that KIC 10920273 and 
KIC 11395018 are at the faint end of the Kepler astero- 
seismic targets; hence, this work sets a lower limit to the 
quality of information we can expect from asteroseismol- 

ogy- 

Wc thank the entire Kepler team, without whom these results 
would not be possible. Funding for this Discovery mission is pro- 



vided by NASA's Science Mission Directorate. We also thank all 
funding councils and agencies that have supported the activities 
of KASC Working Group 1. This article is based on observa- 
tions made with the Nordic Optical Telescope (NOT) operated 
on the island of La Palma in the Spanish Obscrvatorio del Roquc 
de los Muchachos. GD gratefully acknowledges financial support 
from the following institutions: Advanced Study Program (ASP) 
of the National Center for Atmospheric Research (NCAR), NASA 
under Grant No. NNX11AE04G (together with MP and TSM), 
The Danish Council for Independent Research; and thanks Y. 
Elsworth, S. Hekker, M. Stgslicki, J. C. Suarez, M. J. Thomp- 
son, and the anonymous referee for useful comments. NCAR is 
partially supported by the National Science Foundation. Fund- 
ing for the Stellar Astrophysics Centre is provided by The Dan- 
ish National Research Foundation. The research is supported 
in part by the ASTERISK project (ASTERoseismic Investiga- 
tions with SONG and Kepler) funded by the European Research 
Council (Grant agreement no.: 267864). AOT acknowledges sup- 
port from Sondcrforschungsbercich SFB 881 "The Milky Way Sys- 
tem" (subproject A5) of the German Research Foundation (DFG). 
1MB is supported by the grant SFRH / BD / 41213 /2007 from 
FCT /MCTES, Portugal. 1MB and MJPFG were supported in 
part by grant PTDC/CTE-AST/098754/2008 from FCT-Portugal 
and FEDER. JM-Z acknowledges the Polish Minstry grant num- 
ber NN203 405139. KB acknowledges the funding support from 
the INAF Postdoctoral fellowship. This research was carried out 
while OLC was a Henri Poincare Fellow at the Observatoire de 
la Cote d'Azur. The Henri Poincare Fellowship is funded by the 
Conscil General des Alpes-Maritimes and the Observatoire de la 
Cote d'Azur. RAG has received funding from the European Com- 
munity's Seventh Framework Programme (FP7/2007-2013) under 
grant agreement No. 269194 (IRSES/ASK). SGS acknowledges 
the support from the Fundagao para a Ciencia e Tecnologia (grant 
ref. SFRH/BPD/47611/2008) and the European Research Coun- 
cil (grant ref. ERC-2009-StG-239953). TLC acknowledges finan- 
cial support from project PTDC/CTE-AST/098754/2008 funded 
by FCT/MCTES, Portugal. WJC acknowledges financial support 
from the UK Science and Technology Facilities Council. We ac- 
knowledge the KITP staff at UCSB for their warm hospitality dur- 
ing the research program "Asteroseismology in the Space Age". 
This research was supported in part by the National Science Foun- 
dation under Grant No. PHY05-51164. 



REFERENCES 



Adelberger, E. G., Austin, S. M., Bahcall, J. N., et al. 1998, 

Reviews of Modern Physics, 70, 1265 
Aerts, C, Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, 

Asteroseismology, Springer, Heidelberg (Springer, Heidelberg) 
Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41 
Alexander, D. R. & Ferguson, J. W. 1994, ApJ, 437, 879 
Angulo, C, Arnould, M., Rayet, M., et al. 1999, Nuclear Physics 

A, 656, 3 

Appourchaux, T., Chaplin, W. J., Garcia, R. A., et al. 2012, 
A&A, 543, A54 

Bahcall, J. N. & Pinsonneault, M. H. 1992, Reviews of Modern 
Physics, 64, 885 

Bahcall, J. N., Pinsonneault, M. H., & Wasserburg, G. J. 1995, 

Reviews of Modern Physics, 67, 781 
Barnes, S. A. 2003, ApJ, 586, 464 
Barnes, S. A. 2007, ApJ, 669, 1167 

Bedding, T. R. 2012, in Astronomical Society of the Pacific 
Conference Scries, Vol. 462, Progress in Solar/Stellar Physics 
with Helio- and Asteroseismology, ed. H. Shibahashi, 
M. Takata, & A. E. Lynas-Gray, 195 

Benomar, O., Bedding, T. R., Stello, D., et al. 2012, ApJ, 745, 
L33 

Bodenheimer, P. 1965, ApJ, 142, 451 

Bohm-Vitense, E. 1958, Zeitschrift fur Astrophysik, 46, 108 
Borucki, W. J., Koch, D., Basri, G., ct al. 2010, Science, 327, 977 
Brandao, I. M., Dogan, G., Christensen-Dalsgaard, J., et al. 2011, 
A&A, 527, A37 

Campante, T. L., Handberg, R., Mathur, S., ct al. 2011, A&A, 
534, A6 

Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550 
Chaplin, W. J., Appourchaux, T., Elsworth, Y., et al. 2010, ApJ, 
713, L169 



Chaplin, W. J., Elsworth, Y., Isaak, G. R., Miller, B. A., & New, 

R. 1999, MNRAS, 308, 424 
Chaplin, W. J., Kjeldsen, H., Christensen-Dalsgaard, J., et al. 

2011, Science, 332, 213 
Christensen-Dalsgaard, J. 2008a, Ap&SS, 316, 13 
Christensen-Dalsgaard, J. 2008b, Ap&SS, 316, 113 
Christensen-Dalsgaard, J., Dappcn, W., Ajukov, S. V., et al. 

1996, Science, 272, 1286 
Creevey, O. L., Dogan, G., Frasca, A., et al. 2012, A&A, 537, 

Alll 

Deheuvels, S., Garcia, R. A., Chaplin, W. J., et al. 2012, ApJ, 
756, 19 

Deheuvels, S. & Michel, E. 2011, A&A, 535, A91 
Eggenberger, P., Meynet, G., Macdcr, A., et al. 2008, Ap&SS, 
316, 43 

Eggenberger, P., Meynet, G., Macdcr, A., ct al. 2010, A&A, 519, 
A116 

Epstein, C. R. & Pinsonneault, M. H. 2012, arXiv: 1203. 1618; 

submitted to ApJ 
Ferguson, J. W., Alexander, D. R., Allard, F., ct al. 2005, ApJ, 

623, 585 

Fernandes, J. & Monteiro, M. J. P. F. G. 2003, A&A, 399, 243 
Frandscn, S. & Lindberg, B. 1999, in Astrophysics with the NOT, 

ed. H. Karttunen & V. Piirola, 71 
Gelly, B., Lazrek, M., Grec, G., et al. 2002, A&A, 394, 285 
Gilliland, R. L., Jenkins, J. M., Borucki, W. J., et al. 2010, ApJ, 

713, L160 

Grevesse, N. & Noels, A. 1993, in Origin and Evolution of the 

Elements, ed. S. Kubono & T. Kajino, 14 
Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 
Kjeldsen, H., Bedding, T. R., & Christensen-Dalsgaard, J. 2008, 

ApJ, 683, L175 

Koch, D. C, Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 



10 



Dogan et al. 



Krishnamurthi, A., Pinsonncault, M. H., Barnes, S., & Sofia, S. 

1997, ApJ, 480, 303 
Lazrek, M., Baudin, F., Bertello, L., et al. 1997, Sol. Phys., 175, 

227 

Lebreton, Y., Auvergne, M., Morel, P. J., & Baglin, A. 1993, in 
Astronomical Society of the Pacific Conference Series, Vol. 40, 
IAU Colloq. 137: Inside the Stars, ed. W. W. Weiss & 
A. Baglin, 474 

Mamajek, E. E. & Hillenbrand, L. A. 2008, ApJ, 687, 1264 
Mathur, S., Handberg, R., Campante, T. L., et al. 2011, ApJ, 
733, 95 

Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 
152 

Meibom, S., Mathieu, R. D., Stassun, K. G., Liebesny, P., & 

Saar, S. H. 2011, ApJ, 733, 115 
Metcalfe, T. S., Chaplin, W. J., Appourchaux, T., et al. 2012, 

ApJ, 748, L10 

Metcalfe, T. S., Creevey, O. L., & Christensen-Dalsgaard, J. 

2009, ApJ, 699, 373 

Metcalfe, T. S., Monteiro, M. J. P. F. C, Thompson, M. J., et al. 

2010, ApJ, 723, 1583 

Michaud, G. & Proffitt, C. R. 1993, in Astronomical Society of 
the Pacific Conference Series, Vol. 40, IAU Colloq. 137: Inside 
the Stars, ed. W. W. Weiss & A. Baglin, 246-259 

Miglio, A., Montalban, J., & Maceroni, C. 2007, MNRAS, 377, 
373 



Morel, P. 1997, A&AS, 124, 597 
Osaki, J. 1975, PASJ, 27, 237 

Parker, P. D. M. 1986, in Physics of the Sun. Volume 1, ed. 

P. A. Sturrock, T. E. Holzer, D. M. Mihalas, & R. K. Ulrich, 

Vol. 1, 15-32 
Pinsonneault, M. 1997, ARA&A, 35, 557 

Pinsonneault, M. H., An, D., Molenda-Zakowicz, J., et al. 2012, 
ApJS, 199, 30 

Pinsonncault, M. H., Kawaler, S. D., & Demarque, P. 1990, 

ApJS, 74, 501 
Proffitt, C. R. & Michaud, G. 1991, ApJ, 380, 238 
Randich, S. 2010, in IAU Symposium, Vol. 268, IAU Symposium, 

ed. C. Charbonnel, M. Tosi, F. Primas, & C. Chiappini, 

275-283 

Rogers, F. J. & Nayfonov, A. 2002, ApJ, 576, 1064 
Scuflaire, R., Montalban, J., Thcado, S., ct al. 2008, Ap&SS, 316, 
149 

Sestito, P. & Randich, S. 2005, A&A, 442, 615 
Skumanich, A. 1972, ApJ, 171, 565 

Suarez, J. C, Goupil, M. J., Reese, D. R., ct al. 2010, ApJ, 721, 
537 

Woitaszek, M., Metcalfe, T., & Shorrock, I. 2009, in Proceedings 

of the 5th Grid Computing Environments Workshop, p. 1-7 
Zappala, R. R. 1972, ApJ, 172, 57 



TABLE 2 

Input physics used in the evolution codes 



Team 


Diffusion 


Convection 


Overshoot 


EOS 


Opacities 


Nuclear reaction 




& settling 


treatment 


(core) 




(high/low temperature) 


rates 


AMP a (ASTEC) b 


He c 


MLT d 


no 


OPAL2005 C 


OPAL' /Alexander & Ferguson ( 1994) 


B & P (1992)8 


ASTEC1 


none 


MLT 


no 


OPAL2005 


OPAL /Ferguson et al. ("2005) 


NACRE 11 


ASTEC2 


He & heavy elements 


MLT 


no 


OPAL2005 


OPAL /Ferguson et al. (2005) 


NACRE 


CESAM 1 


none 


CGJVP 


yes 


OPAL2005 


OPAL /Alexander & Ferguson f 19941 


NACRE 


Geneva 


He & heavy elements 1 


MLT 


yes 


OPAL2005 


OPAL/ Alexander & Ferguson (1994) 


NACRE 



a) Metcalfe ct al. (2009), b) Aarhu s St ellar Evolution Code { Christcnscn-Dalsgaard 2008a), c) as des cribed by Q vlichaud & Proffitt (1993), d) Mixing length theory dBohm-V itcnsc 
1958), c) Rogers & Navfonov (2002), f) Iglcsias & Rogers (1996), g) Bahcall & Pinsonncault (1992), h)Angulo ct al. (1999), i) Morel (1997), j) Canuto-Goldman-Mazzitclli model for 
turbulent convection l lCanuto et al.111999) . k) lEggenberger ct al. (2008) 1) Proffitt & Michaud ( 1991) 



TABLE 3 

Parameter space searched by each team 



Team 


M/M fi 


Z/X 


Y 




a 




a ov 


AMP (ASTEC) 


0.75-1.75 


0.0026-0.079 


0.22- 


-0.32 


«MLT= 


=1.0-3.0 


N/A 


ASTEC1 


1.00-1.60 


0.01-0.07 


0.24- 


-0.32 


"MLT= 


=1.8 


N/A 


ASTEC2 


1.2-1.4 


0.025-0.046 


0.26- 


-0.30 


"MLT= 


=1.78-1.84 


N/A 


CESAM 


N/A* 


0.026-0.042 


0.24- 


-0.28 


"CGM 


== 0.52-0.68 


0.0-0.2 


Geneva 


1.00-1.50 


0.016-0.040 


0.25- 


-0.30 


«MLT = 


=1.8 


0.1 



* For each given set of parameters (Z/X, Y, a, a v), the method proposed by D eheuvels . 
separation and the frequency of the mixed modes. 



Michel ( 2011) results in a precise estimate of the mass by using the observed large frequency 
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TABLE 4 

Fitted parameters for KIC 10920273. 



Model 


M/M 


(Z/X)i 


Yj 


a 


i(Gyr) 


L/Lq 


R/Rq 


T cff (K) 


log g 


[Fc/H] 


x c 


x 2 ■ 


Aatm 


SA1 (AMP) 


1.00 


0.0154 


0.296 


2.04 


6.74 


3.31 


1.779 


5844 


3.937 


—0.203 


0.0 


6.33 


2.03 


SA2 (AMP) 


1.03 


0.0143 


0.267 


1.98 


7.28 


3.24 


1.797 


5781 


3.942 


-0.235 


0.0 


6.35 


2.10 


SA3 (AMP)" 


0.96 


0.0152 


0.311 


1.94 


6.96 


3.14 


1.754 


5805 


3.932 


-0.208 


0.0 


3.09 


1.90 


SA4 (AMP) b 


1.00 


0.0151 


0.285 


2.02 


6.86 


3.22 


1.780 


5799 


3.937 


-0.210 


0.0 


3.60 


1.85 


SB1 (ASTEC1) 


1.02 


0.0200 


0.290 


1.80 


7.59 


2.97 


1.788 


5674 


3.944 


-0.088 


0.0 


11.44 


1.71 


SB2 (ASTECl) b 


1.06 


0.0300 


0.300 


1.80 


7.45 


2.84 


1.813 


5572 


3.948 


0.088 


0.0 


23.60 


4.21 


SC (ASTEC2) 


1.10 


0.0214 


0.285 


1.81 


5.76 


3.42 


1.824 


5819 


3.960 


-0.090 


0.0 


33.90 


0.79 


SD (CESAM) b 


0.97 


0.0152 


0.300 


0.56° 


7.80 


3.10 


1.760 


5779 


3.930 


-0.207 


0.0 


4.00 


1.90 


SE (Geneva) 


1.15 


0.0220 


0.275 


1.80 


5.60 


3.77 


1.861 


5900 


3.959 


-0.040 


0.0 


45.78 


1.40 


weighted mean 


1.00 


0.0161 


0.294 




7.12 


3.18 


1.776 


5787 


3.937 


-0.188 








standard deviation 


(0.04) 


(0.0030) 


(0.014) 




(0.47) 


(0.13) 


(0.021) 


(55) 


(0.007) 


(0.064) 








Creevey et al. 


1.25±0.13 








4.5±1.8 


3.6±1.2 


1.90±0.05 




3.97±0.03 










(2012) 










5.0±1.9 d 








3.94±0.03° 











a Maximal-list frequencies are used as input. 

b Two frequencies are excluded from the maximal-list frequencies (see text). 

c Canuto-Goldman-Mazzitelli (CGM) model for turbulent convection (Canuto ct al. 1996) is used in this model. 
d This is the value when the average small frequency separation is also used as a sei smic constraint. 

6 Asteroseismic logg obtained from scaling relations (Table 7 of lCreevev et aLll20T2l) ^ 

O 
crq 

! 

TABLE 5 



Fitted parameters for KIC 11395018. 



Model 


M/Mq 


(Z/X)i 


Yi 


ot 


t(Gyr) 


L/L 


R/R e 


T ctf (K) 


log g 


[Fc/H] 


x c 


2 


v 2 

Aatm 


BA1 (AMP) 


1.23 


0.034 


0.301 


1.94 


4.46 


4.40 


2.158 


5697 


3.860 


0.144 


0.0 


8.07 


0.49 


BA2 (AMP) 


1.32 


0.028 


0.241 


1.94 


4.84 


4.53 


2.210 


5671 


3.869 


0.049 


0.0 


8.51 


0.69 


BA3 (AMP)° 


1.26 


0.034 


0.294 


2.02 


4.28 


4.64 


2.175 


5749 


3.863 


0.140 


0.0 


7.36 


0.55 


BB1 (ASTEC1) 


1.235 


0.033 


0.282 


1.80 


5.05 


4.04 


2.156 


5573 


3.861 


0.114 


0.0 


17.03 


1.02 


BB2 (ASTECl) a 


1.22 


0.040 


0.310 


1.80 


4.65 


4.11 


2.154 


5605 


3.858 


0.213 


0.0 


17.08 


1.02 


BC1 (ASTEC2) 


1.25 


0.034 


0.297 


1.81 


4.29 


4.42 


2.265 


5615 


3.820 


0.116 


0.0 


47.34 


0.90 


BC2 (ASTEC2) b 


1.32 


0.043 


0.270 


1.81 


4.89 


4.06 


2.211 


5515 


3.870 


0.243 


0.0 


17.65 


2.01 


BD (CESAM) 


1.29 


0.030 


0.260 


0.64° 


4.50 


4.89 


2.190 


5806 


3.866 


0.026 


0.0 


4.04 


1.19 


BE (Geneva) 


1.32 


0.033 


0.275 


1.80 


4.30 


5.18 


2.207 


5867 


3.871 


0.120 


0.0 


115.69 


1.37 


weighted mean 


1.27 


0.033 


0.276 




4.57 


4.54 


2.184 


5706 


3.863 


0.103 








standard dev. 


(0.04) 


(0.004) 


(0.022) 




(0.23) 


(0.30) 


(0.024) 


(92) 


(0.008) 


(0.070) 








Creevey et al. 


1.37±0.11 








3.9±1.4 


4.2±1.1 


2.23±0.04 




3.88±0.02 










(2012) 










4.5±0.5 d 








3.86±0.03 c 











a Maximal-list frequencies are used as input. 

b No diffusion is included in this model unlike other ASTEC2 models. 

c Canuto-Goldman-Mazzitelli (CGM) model for turbulent convection (Canuto ct al. 1996) is used in this model. 
d This is the value when the average small frequency separation is also used as a sei smic constraint. 
c Asteroseismic logg obtained from scaling relations (Table 7 of Creevey ct al. 2012) 



